import numpy as np
import matplotlib.pyplot as plt

x = np.linspace(-1.8, 1.8, 10000)
y = x**(2/3) + 0.75 * np.sqrt(3.3 - x**2) * np.sin(4.47 * np.pi * x)

plt.plot(x, y, color='red')
plt.fill_between(x, y, color='red', alpha=0.3)
plt.title('Heart Shape Function')
plt.show()